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Quantum entanglement can manifest itself in the narrowing of wavepackets. We define the phe- 
nomenon of phase entanglement and describe its effect on the interpretation of spatial localization 
experiments. 
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Substantial progress has been made in the last decade 
in the understanding of quantum entanglement in contin- 
uous spaces. For example, criteria are now available that 
determine the existence of entanglement, or not, of bipar- 
tite states P, 0, • It is only in continuous spaces that 
the degree of entanglement is not bounded, and it is here 
that the interest in high entanglement is focused. This is 
where the quantum-classical border domain || is entered 
by the Einstein-Podolsky-Rosen (EPR) thought experi- 
ment and examination of high entanglement realms 
associated with EPR physics has begun Q. A substan- 
tial theoretical clarification was made earlier by Reid and 
Drummond || in the context of two-mode photonic field 
amplitudes, which can formally play the roles of posi- 
tion and momentum variables. More recent theoretical 
studies have explored various physical realizations anal- 
ogous to the EPR scenario |3, llfj, HH ■ Experimental 
advances include early observations of small degrees of 
continuum entanglement [l^l IbH fl5j | . Recent measure- 
ments have been made with down-conversion photons in 
the high entanglement domain ^(|, where the EPR un- 
certainty relation deviates strongly from the Heisenberg 
relation, and massive-particle observations are being pro- 
posed [l7j . 

While momentum and position offer equivalent bases 
in the usual mathematical Fourier sense for discussions 
of quantum entanglement in continuous spaces, they fre- 
quently correspond to experimentally quite different op- 
portunities or challenges. In this Letter we show that the 
operational effect of this difference is to expose an unex- 
pected feature that we believe has not previously been 
noted. This feature is phase entanglement. We define 
it and show that it should be observable in a massive- 
particle EPR experiment. 

Surprisingly, one can find phase entanglement even 
in Gaussian states |isj that describe very uncompli- 
cated physical situations. Perhaps equally surprising, 
the original EPR scenario of fragmentation or "breakup" 
of two massive particles is one of these. Let us con- 
sider a breakup process in which the two particles, ini- 
tially joined as a single unit, are decomposed into non- 
interacting fragments that fly apart. Molecular dissoci- 
ation provides a natural physical example. We can an- 
alyze the dynamics of this simple two body problem by 



using center of mass and relative coordinates, in terms 
of which the joint wave function at breakup is separa- 
ble: ^(xi,X 2 ;t = 0+) = <pcm(Xl + X 2 ) <f>rel(xi ~ X 2 ), 

where we have taken equal masses for simplicity in defin- 
ing coordinates. However, there is no guarantee that the 
separability extends to the individual particles, and gen- 
erally they are in fact entangled and ^(xi, x 2 ; t) cannot 
be written in the product form $i(xi,t) & 2 (x 2 ,t). 

To be specific, suppose that the joint state initially 
takes the following form in position space 
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where a > and b > 0. We will consistently omit ir- 
relevant normalization factors. The first factor, the cm 
part of the wavefunction, is a free-particle Gaussian. The 
choice of the exact form of the relative part only affects 
quantitative details but not the interesting qualitative 
features (see 0D> so we have taken it to be Gaussian 
as well. Note that the entanglement in Eq. (|TJ is con- 
trolled by the packet width parameters a and b. The two 
particles are initially uncorrelated (the joint wave func- 
tion factors) when the cross term x±x 2 is absent from the 
exponent, which (for Gaussians) is automatically true in 
the special case when a = b, so we assume a =/= b. The 
product of two Gaussian functions is also Gaussian, of 
course, and many analyses of Gaussian states are avail- 
able, and a great deal is known about conditions on their 
entanglement [Til flij . I20L l2l| . We are raising a new point 
here. 

After breakup the two particles are massive and free, so 
their positional wavepackets spread when they separate 
from each other. The nature of their free space evolution 
is well known. It will be useful to reproduce it here, 
approaching it first from the Fourier transform of Eq. lfl|l: 
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where pi = hki is the momentum of particle i. It evolves 
with the free space Hamiltonian: 
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Since the momenta p\ and p2 are quantum non- 
demolition (QND) operators for the free space Hamil- 
tonian, the time evolution only adds a phase to 'J. We 
have 



where Q is the quantum diffusion length: 
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By Fourier transform we quickly obtain the correspond- 
ing coordinate-space result: 



j_ a 2 (x 1 + 3 ; 2 ) 2 +t, 2 (x 1 - 3 i 2 ) 2 +2iq(xj+x^) 
'4 a 2 b a_ Q 4 + iQ(a a +b a ) 
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We re-emphasize that these Gaussians are not those for 
harmonic oscillators but rather for free particles. The 
time dependence of the ever-growing quantum diffusion 
length Q is important. 

In order to comment on the quantum information con- 
tained in this time-dependent breaku p st ate we will in- 
troduce the Schmidt decompositions [22] of ^(ki,k 2 ;t) 
and ^(xi,X2',t). Since Eq. is similar to a two- mode 
squeezed state, the Schmidt basis here is also spanned by 
number states. Indeed, using identities for Hcrmite func- 
tions, it is easy to find that the Schmidt modes r) 
are given by 



Mk, t) = -W -L H n (K)e-^e-i^, (7) 
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and they determine the characteristic single-sum Schmidt 
expansion: 
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The number of terms in the Schmidt sum is infinite, 
but since the A ra 's are reduced density matrix eigenval- 
ues we have Y) X n = 1, and there is generally a small 
number that contribute significantly. This number is well 
estimated by what is reasonably called the Schmidt num- 
ber K, given by 
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where we defined e r = a/6. Similar expressions have been 
found useful in a study of transverse high entanglement 
in down conversion 101. As we could have expected, 



^> has the same Schmidt form as a two-mode squeezed 
photon state, and the squeezing parameter is r. One has 
1 < if < oo, so if is a suitable choice for a quantitative 
measure of entanglement, and is also experimentally of 
direct relevance since it defines how many independent 
modes, i.e., Schmidt functions (f> n (K,r), are important in 
any given study. 

Evaluations of entanglement can be associated with 
variances of Gaussian states in several ways 0, H S HH • 
Actual measurements of the degree of continuous entan- 
glement are just beginning to be carried out in the do- 
main of very high entanglement where one approaches 
the EPR limit of correlation. A high-entanglement mea- 
surement scenario for a pair of particles undergoing an 
EPR-type breakup has been described for the pho ton and 
recoiling atom in stimulated light scattering [ll| and for 
the signal and idler photons in parametric down conver- 
sion [l0( for which the first related experimental results 
have recently appeared [r| ■ A theoretical proposal for a 
similar experimental determination, but for massive par- 
ticles, was introduced in Ref. 01 > an d a cold optical lat- 
tice context for such an experiment has also been identi- 
fied ^3 • These situations are of wide interest because in 
the massive-particle case one will reproduce for the first 
time the exact EPR scenario. 

In all of these studies of EPR-like entanglement the 
roles played by conditional and unconditional variances 
are central || . In momentum space, the calculations are 
straightforward for conditional and unconditional vari- 
ances for either of the two particles. They correspond 
to measurements of dispersion in momentum that are ei- 
ther single-particle or coincidence measurements. Since 
the momentum space time dependence appears only in 
the phase of the two-particle wave function Q, the mo- 
mentum variances do not change as time evolves. We 
find for the single-particle dispersion 
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On the other hand, a coincidence measurement gives the 
conditional dispersion: 
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where the expectation values are now calculated with 
respect to the conditional probability P(fci|fc 2 ) = 
P(fci , k?) I P(k?.). By combining these in the Fedorov ra- 
tio [12j we obtain 
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This expression shows that the ratio of single to coin- 
cidence packet widths R p , which is experimentally di- 
rectly accessible, is a quantitative entanglement measure 
through its connection with the Schmidt number if. 
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Under some experimental conditions, it is easier to 
measure position variances. We evaluate them for our 
wave function and find: 
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where the expectation (• • • ) X2 is with respect to the po- 
sition conditional probability P(xi|x2), and we define 
to = mab/h. 

Entanglement is representation-independent (and 
time-independent) in such an interaction-free context as 
we are considering. Thus, given the simple Fourier- 
inverse character of momentum and position representa- 
tions, it is surprising that R x = Axf nslc (t)/Ax c 1 oinc (x 2 ;t) 
does not exhibit the same properties as R p . Not only 
that, the x-space width ratio is time-dependent: R x — 
K C(t), where 
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A plot of the function C(t) is shown in Fig. ^ 

The unexpected distinction between R p and R x reveals 
a novel element of free-particle Gaussian states (in con- 
trast to the more commonly studied harmonic oscillator 
Gaussian states) that is open to experimental study, as 
we now explain. If we inspect the probability density in 
position space (ignoring irrelevant normalization factors): 



P(x u x 2 ;t) = \^(x u x 2 ;t)\' 
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we easily see that only the first factor carries the en- 
tangling x\x 2 product term in the exponent. Therefore 
evidence of entanglement vanishes when the first expo- 
nent vanishes, i.e., when either (i) a = 6, or (ii) Q 2 = ab. 
The first case is trivial as remarked earlier. The second 
provides the new feature of Gaussian entanglement. It 
identifies the time t — to at which entanglement informa- 
tion in ^(xi,X2',t) in Eq. |jfj| is completely transferred 
to the phase of the state. At such a time none of the 



1.2 
1 

0.8 

C(0 0.6 
0.4 
0.2 




r = 










\r = 2 






= 5 





10~ 3 10' 2 10" 1 1 10 10 2 
t/t 



w 3 



FIG. 1: Plot of the function C{t) in Eq. TOl at r = 0. 1. 2 
and 5. The time is in units of to. 



coherences in the two particle density matrix appear in 
the joint probability density 11711 . 

Moreover, this result makes it clear that not only at 
t = to but at all times t ^ to some portion of the en- 
tanglement is lost from the spatial joint probability, i.e., 
has become entanglement in positional phase alone. The 
function C (t) is a natural measure of the extent to which 
entanglement information has been shifted into the phase 
of the two-particle state. Clearly we have the ordered in- 
equality: 
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and C(t) shows that only two times, t = and t — > oo, 
provide R x = R p = K . The evolution of R x away from 
and then back to the value R p occurs for an interesting 
reason. It is the means by which the quantum state en- 
sures a smooth transfer between two physically different 
conditions of breakup entanglement. 

These two conditions are easiest to see as distinct in 
the case of high entanglement, which occurs if a ^> b 
or b^$> a. In either case the packet that initially has the 
much smaller width (cither center of mass or relative) has 
to spread much faster than the other packet (Heisenberg 
uncertainty principle in action) . Inevitably, what was the 
wider but slower-spreading packet becomes the narrower 
(Einstein derealization and relocalization in action) and 
their roles in describing entanglement become reversed. 
Asymptotically their final ratio is just the reciprocal of 
their initial ratio, but as our expression 1)10(1 indicates, a 
reversal of the width ratio simply returns K to the same 
value it had initially. We also can see that the trivial 
pure-phase condition a = b also fits, because when the 
widths are initially equal they spread at the same rate 
and no Einstein localization ever takes place. 

The Einstein conditional uncertainty and the uncondi- 
tional Heisenberg uncertainty can now be compared for 
a general breakup scenario. The latter uncertainty is of 
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course preserved throughout the coordinate-space shape- 
changing that occurs: 
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for all t and r. On the other hand, the conditional Ein- 
stein uncertainty relation reads 
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In summary we have examined both momentum and 
position variances arising in a breakup experiment in- 
volving two massive particles, in the case that the par- 
ticles exert little or no influence on each other following 
the breakup. This is a situation that can be realized 
in several ways, including molecular dissociation |l2| . 
Momentum and position are formally equivalent bases 
from which to examine entanglement, while they corre- 
spond to experimentally different scenarios. We showed 
that they also have a fundamental difference in phase en- 
tanglement which becomes evident in comparing the Fe- 
dorov width ratios for momentum and coordinate vari- 
ances. This can arise only in the massive free-particle 
coordinate-space case, exactly the situation of the origi- 
nal Einstein-Podolsky- Rosen thought experiment. Phase 
entanglement has not been discussed previously. We pre- 
dict that it will be detectable in experiments designed to 
achieve high values of entanglement, of fundamental in- 
terest because of their nearness to the classical-quantum 
border [5j- 

Finally, it can be noted that there exist pure phase 
entangled states, e.g., 



~ exp {—/i 2 (x\ + x 2 ) + iv 2 x\X2\ 



(21) 



This is very similar to an EPR state, except that X\ is 
tightly related to k 2 instead of x 2 , as can be seen by 
taking a partial Fourier transform of ^ with respect to 
X2 ■ For such a state the two-particle probability shows no 
correlation between the two particles, and entanglement 
tests based only on variances of position or momentum 
are unable to reveal its entangled nature. 
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